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Abstract 

A system of nonlocal electron-transport equations for small perturbations in a magnetized plasma 
is derived using the systematic closure procedure of V. Yu. Bychenkov et al, Phys. Rev. Lett. 
75, 4405 (1995). Solution to the linearized kinetic equation with a Landau collision operator is 
obtained in the diffusive approximation. The Fourier components of the longitudinal, oblique, and 
transversal electron fluxes are found in an explicit form for quasistatic conditions in terms of the 
generalized forces: the gradients of density and temperature, and the electric field. The full set of 
nonlocal transport coefficients is given and discussed. Nonlocality of transport enhances electron 
fluxes across magnetic field above the values given by strongly collisional local theory. Dispersion 
and damping of magnetohydrodynamic waves in weakly collisional plasmas is discussed. Nonlocal 
transport theory is applied to the problem of temperature relaxation across the magnetic field in 
a laser hot spot. 

PACS numbers: 52.25.Fi, 52.35.Bj, 52.38.Fz 
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I. INTRODUCTION 



Nonlocal transport theories remain the subject of active studies in plasma physics espe- 
cially in laser produced plasmas. The strong dependence of the coUisional mean-free-path 
on the particle kinetic energy limits the validity of standard local transport relations to very 
long inhomogeneity scale lengths L. The value of L has to be at least a hundred times 



longer t 
proach 



lan the electron-ion mean- free-path, Agj, for the Chapman-Enskog perturbative ap- 
to be valid P|. This condition can be easily violated in most of the hot plasmas 



produced by high energy lasers. Likewise, many magnetized plasmas in laboratory ^ and 
astrophysics require kinetic modeling because of the large values of \eilL. In the presence 
of a strong magnetic field electron transport can be localized in the transverse direction by 
the short gyroradius, p. Therefore, often only nonlocal transport in response to longitudinal 
gradients is considered (cf. e.g. j^). We will show that for an arbitrary direction of the 
gradient, longitudinal transport coefficients are affected by the magnetic field in the nonlocal 
transport regime. 

Nonlocality of thermal transport has been well documented in experiments. For example, 
no agreement was found between theory and experiment for time-resolved two-dimensional 
images of x-ray emissions 0, 0] until nonlocal electron transport was included. It is also clear 
that self-generated magnetic fields play an important role in the evolution of laser fusion 
plasmas ITI. These observations have motivated our current study, which also complements 
previous works on magnetized plasmas ^ a, in the regime of weak collisionality. 

We have applied the method of Ref. 10| to derive nonlocal closure relations for a plasma 
with a magnetic field. The result is equivalent to the linearized kinetic theory of small 
amplitude perturbations. Our derivation involves a solution to the initial value problem 
for the electron kinetic equation with electron-electron (e-e) and electron-ion (e-i) Landau 
collision operators. The linearized kinetic equation is solved in the Fourier k-vector space 
for slowly varying processes. We use approximations which make our results correct for 
high-Z plasmas. We have neglected e-e collisions in the equation for the anisotropic part of 
the electron distribution function. Electron-electron collisions are only kept in the equation 
for the isotropic part of the distribution function which is unaffected by e-i collisions. In 
addition, the anisotropic part of the electron distribution function (EDF) is approximated 
by a single angular harmonic. This corresponds to the so called linearized diffusive Fokker- 
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Planck model 
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12l lisf. The original study 13] of an unmagnetized plasma employed 
a summation technique [isl. Q], which included an infinite set of angular harmonics and 
thus ensured validity of the transport theory over the entire range of particle collisionality 
including collisionless limit. Here, because of the geometrical complexity introduced by the 
magnetic field we are unable to carry the summation over the infinite set of harmonics. 
Therefore our results are applicable to an inhomogeneity scale length corresponding to k- 
vectors such that fc_Lp < 1 and k\\Xei < 1, where k± and fcy are wave vector components with 
respect to the direction of the magnetic field. However, in the case where the spatial gradients 
are almost orthogonal to the direction of magnetic field our theory is valid even for kXgi > 1 
provided kX^i < Xei/p- The validity condition A;||Aei < 1 is more relevant to dense plasmas 
as for example in inertial confinement fusion targets. Finally the closure relations derived 
in this paper apply to small amplitude perturbations of temperature, density and electric 
field about the homogeneous equilibrium state. A rigorous generalization of these results 
to nonlocal and nonlinear transport theory is a difficult problem. One possible solution, 
Ql is based on linear theory results for varying in space background parameters. One can 
also gain a valuable insight into validity of different phenomenological models reproducing 
Fokker-Planck simulation results as in Ref. 

An important difference between the nonlocal theory of unmagnetized plasmas and 
that for magnetized plasmas is that one now deals with three different components of the 
electron fiux: longitudinal, transverse, and oblique. All of these are calculated here in 
the Fourier-space and are related by closure relations to perturbations in the temperature, 
density and electric field. Hydrodynamical variables have been defined in such a way that in 
strongly collisional limit we recover classical local results of the Chapman-Enskog method 

Q. 

The paper is organized as follows. Section ITTl presents a linearized solution to the kinetic 
FP equation for the EDF in a magnetized plasma. In Sec. IIIII we use this solution to 
derive the nonlocal electron fiux. The strongly collisional limit is examined in Sec. IIVI 
Section |3 is devoted to study of nonlocal transport coefficients and their dependence on 
the collisionality parameter and the Hall parameter. Nonlocal hydrodynamic equations for 
small perturbations are formulated in Sec. |VT1 and are used to derive dispersion relations 
and damping of MHD waves. In Sec. I VI II we apply our transport theory to study of the hot 
spot temperature relaxation in the context of laser produced plasma with a self-generated 
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magnetic field. We conclude with a discussion and summary in Sec. IVIIII 



II. ELECTRON DISTRIBUTION FUNCTION IN A MAGNETIZED PLASMA 

We start from the standard electron kinetic equation for the EDF, /e, in a high-Z plasma 
placed in a homogeneous magnetic field B 

^ + V- ^ - — fE+ i[v X BlV ^ = CMe] + C^elfeJe] , (l) 

at or rUe \ c J ov 

where E is the electric field, c is the speed of light, e and me are the electron charge and 
electron mass, Cei[fe] and Cee[fe, fe] are electron-ion and electron-electron collision terms 



Following the previous study [l^, we develop a new nonlocal transport theory of a 
magnetized plasma which describes the plasma response to small-amplitude perturbations. 
The EDF, /e, is written as the sum of Maxwellian function, Fq, (with homogeneous back- 
ground electron density and temperature Tg) and an isotropic perturbation, 5/o, and an 
anisotropic perturbation, 6 fan- fe = Fq + Sfo{v) + Sfaniy)- Small amplitude linear pertur- 
bations are described in terms of Fourier components characterized by the wave vector k. 
The linearized kinetic equation ((T)) gives two coupled equations for the perturbed EDFs, 5/o 
and 6 fan, in k-space 

^ + z < k ■ v/,„ > -^(k ■ u)v^ = ae[5/o, Fo]; (2) 

^ + [nX V]^ + ,k ■ v/.„ - Z < k ■ yrfan > -C4fan] = (3) 

IdFo 



eE9Fo 

-zk ■ v/o H — + I 

nip (7v 



(k-v)(u-v)-l(k.u)t;2 



V dv 



where u is the velocity of an ion flow, Q = eB/m^c is the electron gyrofrequency, brackets 
< ... > denote the averaging over angles defining the orientation of the velocity vector v 
and E = E + uxB/c. We have neglected electron-electron collisions in Eq. (jHl) because 
they play a smaller role compared to electron-ion collisions in the evolution of an anisotropic 
distribution function in a high-Z plasma. 

We introduce a unit vector b = B/i? in the direction of the magnetic field. Three 
components of a vector A can be defined with respect to b in the following manner: Ay = 
b(A ■ b), A/\ = b X A, and A_l = b x [A x b]. Assuming long wavelength perturbations 
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such that k_L • V ^ ma.x{Q, z/gj), ky ■ v ^ Uei and slow time variations d/dt <^Vt the solution 
to Eq. Q for the anisotropic part of the EDF reads (cf. Ref. jjol) 



5 fan 

eFo 



k\\v 



kA • V k I • V 

^l— — + 



1 



+ Ve 



E, 



k|| • v)(u|| ■ v) + (k_L • v)(u_L ■ v) 



-k • u + 



(4) 



3t^e»(^||(k± ■ v) + A;||(u^ • v)) - VL{u\\{y.,^ ■ v) + k\\{\i,s, ■ v)) 

3z/e.fi((k^ ■ V)(UX ■ V) - (k^ ■ v)(Ua ■ V)) + 2^]2((k^ . v) (u^ " v) 



(kA ■ v)(uj_ ■ v)) 



3z/ei(9z/2. + 4fi2) 

where Vei = ^iiZnee^K/m?^v^ is the velocity-dependent electron-ion collision frequency, and 
A is the Coulomb logarithm. Expression has been derived in the approximation, which 
corresponds to the so-called diffusive limit of the Fokker-Planck equation [cf. e.g. ([l^l)]. 
This amounts to keeping only leading terms in the angular expansion of 5 fan as discussed 
in Ref. (jbl). Such approximation may always be justified by invoking a sufficiently large Z 
and magnetic field values for the conditions k±Max{p, Agi} < 1 and fcyAej < 1 to be satisfied. 
By substituting 6 fan © into Eq. Q, we obtain an equation for the isotropic perturbation 
5fo 



dt ' 3 

= Cee[Sfo, Fq 



zk ■ u— + 



ie(E ■ k)Fo 



(5) 



kA ■ [E X k]A;||^] 
k'^i'pi 



k - E, 



To derive nonlocal transport relations, we follow the method of Ref. Q| and first solve an 
initial value problem for Eq. In a weakly collisional plasma the perturbation at t = 
provides a driver for transport processes. We postulate the initial 5/o to be the linearized 
Maxwellian EDF, which is defined by density and temperature perturbations 5ne \t=o= 5n(0) 
and 5Te \t=Q= 5^(0), 

n(0) ^r(O) f 3' 



5/o(t^,0) 



Fo{v) 



(6) 



With the initial condition (jUI), a solution to Eq. for 5/o is expressed in terms of the basic 
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functions ^I/"^: 

Fq k^Te Tie 2 Te 

I ^ek^ • [E X k]A;|| iek • 

where 6fo ((Tj) depends on four parameters: 6n{0), 5T(0), u and E. The ion flow velocity 
can be determined from a simple hydrodynamical model for cold ions and the field E is 
usually evaluated from the quasineutrality condition. The solution given by Eq. (|7j) allows 
elimination of initial perturbations 6n{0), 6T{0) in terms of their current values by taking 
moments of 6fo. To do so and to proceed with closure relations for higher order moments we 
must find basic functions which satisfy the following equation with different source terms, 
Sa: 



m^-^C,Jym^,F,] = SA. (8) 

-TO 



1^1 + ^\ 

Sources on the right hand side of Eq. (jHl) are defined as follows Sn = 1,St = f ^/3f|.g — 1, 
Su = Sn + St, Sei = — S±, = •S'a and Sa = k'^v'^/Si'a, where the subscript a assumes 
the values ||, A, and ±. The effective longitudinal and transverse collision frequencies read 
as follows 

Z/|| = Uei , iy± = + n^)/Uei , Z^A = {^li + ■ (9) 

Expressions Q have been introduced following the standard definitions used in Ref. Q]. 

As already stated, we can use Eq. ((Zj) to calculate hydrodynamic moments such as density 
and temperature perturbations 



5ne = 47r / dvv'6fo, (10) 
Jo 

STe = ^ dvv\v' - Sv'tM ■ 
3ne Jo 

and employ them to eliminate 6n{0) and 6T{0) from Eq. ((Tj). As a result, the EDF is 
expressed in terms of its instantaneous hydrodynamic moments 611^ and 6Tf. and vectors u 
and E as independent generalized thermodynamic forces. Moreover, because of the approxi- 
mation k\\Xei ^ 1 we have found that ion fiow does not contribute to 6fo explicitly. This has 
not been the case for short wavelength perturbations, fcyAej ^ 1 in studies of unmagnetized 
plasmas (cf. ^^). The cancellation of terms containing u in the expression for 5/o occurs 
because of the specific form of equation (|HJ, which defines \E'" as a linear combination of 



and Thus, the distribution function SJq can be expressed in terms of the velocity 
moments of the basic functions, = {An/ue) dvv^"^^ SaFq, as follows: 



Fo 



ie(E ■ k) ^ ie{E* ■ k) Jj^^ - J^^^ 



r)NT 



jjNT 



+ 



ie^^^ ■ [E X kl 



D 



NT „i,N 



r)NT 



^TN iT 

nNT ^ 

^NT 



+ 



iek ■ E/ 

FT 



k^Te 

. 



X 



JJE2T 

^NT I.N 

nNT ^ 

J^NT 



DE2N 
^TN i.T 
r)NT ^ 
^NT 



(11) 



where D^^ = JaJb-JaJr and E* = E + ik(Te/e)f5Te/Te + 5rae/njis the effective electric 



field which has been introduced in both, classical 111 and nonlocal 



id ] transport theories. 



III. NONLOCAL ELECTRON FLUXES 

The expression for an anisotropic part of the EDF ^ is used to calculate electron fluxes: 
the electric current, j, and the heat flux, q, 

] = -ejdwwfe, q=-T,y"dvv(^^-^^/,. (12) 

Terms proportional to the ion velocity, u, in Eq. (0)) do not contribute to electron fluxes 
(I12|) . This is due to the approximation k\\\ei ^ 1. Equations (fT^ ignores in fact small 
contributions oc u x O(fcyAgj). 

Closure relations for the electric current density and the electron heat flux ()12|) have 
the same form as local transport relations in strongly coUisional plasmas. They involve 
linear combination of Fourier components of an effective electric field E* and a temperature 
gradient (VT)k = ikST^, 

j = cr ■ E* + +ia^ ■ k6Te + ■ E , 

q = a'l-W + ik- x5Te + d ■ E . (13) 

However, the two transport matrices, <t and a. have no counterparts in the local relations Q|. 
Both coefficients vanish in the limit of A; — 0. The transport coefficients in (fT^ are: electric 
conductivity tensors cr, thermoelectric tensors o: and temperature conductivity tensors X- 
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They are defined as follows {c = j,q) 




A;2 








\ 



a , 



(Til 



(721 



k\\k I 



-0-12 



0^22 



-ail 



"a . 








-0-21- 



fc2 



k, 



V 



-(Tii- 




( 



OL 



k' 



"11- 



X± -Xa 




(14) 



"21 



fc2 

k\\k I 



-ai2 



"22 



-an 



A;2 

k\\k I 



V 



fc2 



-"21- 
«32-^ «31-^ 



where k = {k±, 0, /cn). Explicit expressions for these coefficients are given in Appendix A. All 
nonlocal transport coefficients ()A2|) are functions of the wave vector and can be calculated 
by using solutions to the equations for the basic functions 

Equations represent the main results of our theory. They define constitutive relations 
for electron fluxes in terms of thermodynamic forces. Transport relations ()13p are linear 
in the temperature gradient and in a generalized field, similarly to local relations of the 
collision dominated hydrodynamics Q]. The nonlocality of expressions fjl3|) is displayed 
in terms of convolutions in the configurational space involving temperature gradient and 
generalized electric field, and kernels which are defined by the set of nonlocal transport 
coefficients. This new form of transport theory generalizes previous results by including 
magnetic field effects. The interplay between weakly collisional and magnetic field effects 
on transport theory is still poorly understood. However it can play an important role in the 
hydrodynamical evolution of laser produced plasmas as was demonstrated by the numerical 
simulations in Ref. j^. 

The transport coefficients ()A2j) depend on the following set of parameters: fcyAej, k±p, 
n/i^^, and Z, where p = vxe/^ is the electron Larmor radius and = A/2/97r Vei{vTe) = 
VTe/^ei is the elcctron-iou collision frequency for thermal electrons In Sec. |3the equiv- 
alent set of independent parameters: kXei (nonlocality parameter), fi/z/^ (Hall parameter), 
k±/k\\, and Z (ionic charge) will be used. 



IV. TRANSPORT THEORY OF A STRONGLY COLLISIONAL PLASMA 

As a benchmark test for our transport theory we first consider the limit of strong collisions, 
i.e. Zk^Xl^ <ti 1. In this limit the basic functions {A = N, T, u, Ei, E2) have the following 
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form 



VTe 



- 1 



(15) 



Expression ()15|) corresponds to a two polynomial approximation to the Laguerre polynomial 
series expansion of the symmetrical part of the EDF P]. Coefficients Cq and are found 
from the solution to Eq. (jSI) in the Lorentz plasma approximation. From the moments of 
the EDF, J^, we can evaluate transport coefficients ()A2jl . 

n 47re^ , o „ n 47re 



Xa 



FT, 



(iff SaFo 



An 



e Jo 



dvv S, 



dvv'^Sg^Fo 



Qa 



24e 



0, 



e Jo 

a = 0, 



(16) 



where Sg^ = k'^v'^iy'^ /v"^^ — b)/Qva {cl =||, A, and 1) and the effective collision frequencies, 
z/q, are given by Eq. (0). With expressions (fTH|) . the electron fluxes ffT^ assume the standard 
form: 



j = ^ CTaEa + ^ aa«ka(5re , 

a a 



(17) 

(18) 



As expected, longitudinal transport coefficients (a =||) do not depend on the magnetic field 
and correspond to the Spitzer-Harm (SH) result without magnetic field. 

In the absence of a magnetic field {Q 0) all transverse and longitudinal transport 
coefficients are equal. The oblique coefficients are proportional to Q, i.e. 



a, 



210enefi 



Xa 



lOlSneTeQ 



(19) 



In the limit of a strong magnetic field Q oo, transverse and oblique transport coefficients 
behave as oc and oc correspondingly. 



0- 



2 T 



a I 



4 



1/31 



t5' 



1/3 



a. 



2mM \ 6fi5 



Xa 



5neTe 



(20) 

(21) 



with e xcep tion of the coefficient which varies as VL~^/'^. Note, that this fractional depen- 
dence |l9| on the electron gyrofrequency follows directly from the integral definition of the 
oblique thermocurrent coefficient in Eq. (fTB|). 







To simplify a comparison with results of classical transport theory jjj, ll9| we rewrite the 
expression for an electron heat flux (fTHj) in terms of the electric current density instead of 
the effective electric field by using the following relation 

E* = -nmj{e^n,)]r. + J] y^Ja - l3fiK5T, . (22) 

a a 

Similarly, one can represent q in the following form 

q = - 5^ Tel3fu - J2 ^MT, . (23) 



The above expressions are now indentical to Ohm's ()22|) and Fourier's ()23|) laws as derived 
in Ref. (the Braginskii's coefficient is replaced here by /Xq). Coefficients /Iq, k^, and 
13^^ in Eqs. ^ and (j2Sl) are defined as follows 

/i|| = -o, /5r = 4' '^ll=X||-A^, (24) 

II II II 

A'-J- ~ n 2 . TP? ' ^ 



(a° a° + «) T,(a<| + 2a0 - «f 

P± - To- — — ' - X± - 



Transport coefficients given by (j24p are the same as the results of Ref. , where a detailed 
comparison with the strongly collisional results of Braginskii [1] was made. Reference 
also gives the functional dependence of and on electron gyrofrequency Q. 



V. NONLOCAL TRANSPORT COEFFICIENTS 



In this section we will describe the dependence of the transport coefficients on the coUi- 
sionality parameter, kXei, for different values of the Hall parameter: f2/z/^ = 0.1 — 100. We 
start with the longitudinal transport coefficients. In the local transport theory the longitu- 
dinal coefficients do not depend on magnetic field. This is not true in general for nonlocal 
transport, except for the special case when the direction of the plasma inhomogeneity is 
along the direction of the magnetic field (A; = k\\). In the general case {k ^ k\\) the depen- 
dence of transport coefficients, including longitudinal coefficients, on the magnetic field is 
defined by the symmetric part of the EDF which is a function of k±p. 
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In the case of weak plasma inhomogeneity along the magnetic field, k± 3> fcy, transport 
coefficients cr|| and x\\ normalized to their classical strongly coUisional values are shown in 
Figs. Hand 121 as functions of the coUisionality parameter kX^i. The nonlocality of transport 
coefficients is characterized by very different functional dependence on kX^i depending on the 
magnitude of the Hall parameter, ^/t^Ji- For small values of the Hall parameter, ^/t^Ji < 1, 
the transport coefficients steadily decrease with increasing kXei- However, as shown in Figs. 
HI and El they can exceed the standard classical values at kXei > 1 for fi/z/^ > 1. Note, that 
the regime of kXei ^ 1 can still be within the validity condition of our theory, k\\Xei < 1, 
when the plasma gradients are almost transverse to the magnetic field, i.e. k± ^ k\\. 

Transport in the case of a strong plasma inhomogeneity along magnetic field k± <^ k\\ is 
similar to the well known zero magnetic field case 10|, therefore the next example considered 
is for k± ~ Figures El and m show transport coefficients a\\ and x\\ functions of kXei 
for k± = k\\. This corresponds to an angle of 45° between the magnetic field vector and the 
direction of a linear gradient. As it is shown in Figs. El and 01 electrical and temperature 
conductivities increase with magnetic field reaching maximum values at f^/t'^ ~ 0.3 (cf. 
curves 3 in both Figs. ElandEj). Further increases in the parameter i^/z/^ results in a small 
decrease in the coefficients cr|| and x\\- At large values of ^/i^J'i, they asymptotically approach 
curves 4 in Figs. OlandlH In general, the longitudinal transport coefficients display smaller 
effects of nonlocality for larger magnetic field values. 

Unlike the longitudinal coefficients, the perpendicular and oblique transport coefficients 
demonstrate strong dependence on magnetic field. To illustrate this, coefficients cta,± and 
Xa,± are shown in Figs. El - IHl as functions of coUisionality parameter for different magnetic 
field values. In general, a strong magnetic field suppresses the effects of nonlocality. For 
instance, the transverse temperature conductivity, X-L? is practically independent of the 
coUisionality parameter over the entire range of variations, kXei < 1, at ^/i^Ji = 3 (curve 
5 in Fig. IH)). As in the longitudinal transport case, the coefficients (Ja,± (Fig. Ej) and 
increase slightly with the coUisionality parameter at l^/z/^?^ ^ 1. One can also see that 
nonlocal effects are stronger for oblique transport coefficients as compared to transversal 
coefficients (cf. Fig. Eland Fig. [71 and Fig. Eland Fig. EJ. 

Figure El shows the ratio of the perpendicular to the longitudinal temperature conductiv- 
ities as a function of Hall parameter. This ratio increases with the value of the nonlocality 
parameter, kX^i- Hence, in a plasma with steep gradients, the inhibition of electron fiuxes 
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due to a magnetic field is not as pronounced as in the strongly collisional case. This corre- 
sponds to an increase in the effective collision frequency (decrease of a mean free path) in 
the nonlocal limit Q]. The kinetic effect of reducing the electron mean free path due to 
transport nonlocality was included in recent 2D hydrodynamical simulations 6] by replacing 
parameter {Q/uJ^"^ with c/(fi/i/^)^ where d is a constant that should be defined from the 
solution to the FP equation. However, our results demonstrate that the phenomenological 
constant d should be more correctly replaced by a function of the gradient scale length, fcAei. 

To complete our discussion of nonlocal, magnetic field dependent transport coefficients, 
note that the analysis of kX^i and fi/z/^ dependencies can be performed for all other coef- 
ficients, a^, 5", and a. These are explicitly calculated in Appendix A. As in the nonlocal 
transport case without a magnetic field 10|, the coefficients in the present theory may 
change signs at some wave numbers. The tensor coefficients, a and a, are negligible as 
f2/z/^ — s> 0, in agreement with the local transport theory |]]]. They are also very small for 
the case of a strong magnetic field ^/i^Ji — > oo. This is the reason why these coefficients 
contribute to the electron fluxes only at fi/z/^ < 1. Finally, all conclusions reached in the 
above discussion are valid only in the limit fcyAej < 1 and k± < Max{l/p, 1/Xei}- These are 
the required conditions for our theory to be valid. 

VI. MAGNETOHYDRODYNAMIC WAVES 

The first two velocity moments of the Eq. (j21) yield continuity and energy balance equa- 
tions for the electrons 

dfi 1 

+ iuek ■ u - -zk ■ j = , (25) 

dn 2 2Te . 2Te., 

— 1 ik ■ q tk ■ J H ik ■ u = (J . 

ot Srig Serif, 3 

Equations ()25|) and the hydrodynamic equation of motion for ions, 

du Ze~ 1 ^ , . 

_E+ R,,, 26 

ot rrii niiTLi 

together with nonlocal transport relations (fTSj) can be used in description of linear per- 
turbations in a magnetized plasma. An ion-electron friction force, Rje, is defined by the 
anisotropic part of the EDF and can be expressed in terms of the electron current and the 
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effective electric field, 

R-je = / dwi^eife = -CngE* - ""'""""' Ja , (27) 



and eventually in terms of the transport coefficients and a. No other transport coefficients 
(c.f. Ref. contribute to Rjg in the limit of kX^i < 1. 

By introducing perturbations of the form oc exp(— icut) into electron ()25p and ion ()26|) 
fluid equations and by using Eq. (jTTj) we obtain 

eujSrie + k ■ j — engk ■ u = 
3e 

— rigCo'^Te — ek ■ q + Tgk ■ j — eTgriek • u = (28) 
uju = kc l + — - t Ja , 

where fij = ZeB/rriiC is the ion gyrofrequency. These equations together with the closure 
relations ()13j) and the standard wave equation, 

k^E- (E-k)k-47ricj/c^j = 0, (29) 

allow a study of linear modes of the magnetized plasma in the low frequency limit, oo <C 
kc, Upe, where Upe is the electron plasma frequency. 

The determinant of the system of equations (j28p and (|29|) gives a dispersion relation for 
quasistatic magnetohydrodynamic (MHD) waves. It has been derived in the Appendix B in 
the following form 



C(cu2 - kc^) + D{uj'^U - clk\) + EVliU = , (30) 



a = J^(^^-A:M), k\ = ^ico^-kivl), 



Symbols used in Eq. ()30p are defined in Appendix B. In the strongly coUisional limit, 
fcAej 0, the dispersion equation (jHUI) has been studied in Ref. '2Cf] in the isothermal 
approximation. In this approximation, long wave contributions given by coefficients A and 
B are neglected because they account for the adiabatic plasma response. Without a magnetic 
field, Eq. (jHIH) corresponds to the ion acoustic dispersion relation (33) of Ref. |2l| . 

Three MHD modes are defined by different branches of the solution to Eq. ()30p for the 
real part of uj. They are shown in Figs. ITUk and II lb and correspond to the well known low 
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frequency waves. The highest frequency branch (curves 1 in Figs. ITUk and ITTk) represents 
the fast MHD wave, u = u+ where ul = {e{v\ + c^) ± ^k^{vl + c2)2 - Ak^kjvlc'i)/2. 
The second branch (curve 2) includes, from small to large values of kXef. the Alfven wave 
{uj = k\\VA <^ ^i), first ion cyclotron wave {uj = Qi) and ion acoustic wave {uj = kcg), 
respectively. Finally, the lowest frequency branch (curve 3) describes the slow MHD wave, 
UJ = uj^ Qi at small kXei and a second ion cyclotron wave {u = Qik\\/k) for large kXei- 

Depending on plasma conditions, nonlocal effects may change the damping of MHD 
modes. This is shown in Figs. ITUb and ITTb for a propagation angle of 80° with respect 
to the direction of a magnetic field. The almost transverse direction of the vector k allows 
examination of the high kX^i values in Figs. and ^2 within the applicability limits of our 
theory, k\\Xei ^ 1. Two parameters, p = {upeVTe/ ^JiC)"^ and h'^ = (Q/uJ^)'^, define different 
regimes of MHD wave damping. The p/h"^ is a well-known /5-parameter, (3 = c^/v^. Figure 
ITUl corresponds to a high-/? (high pressure) plasma and Fig. [TTl describes a low-/? (low 
pressure) case. Note, that we do not discuss effects of ion viscosity. This can be included 
independently by using the standard approach of Ref. [2^. 

In a high-/5 plasma all modes except the fast MHD wave experience the usual collisional 
damping. The damping of the fast MHD wave in this case is unchanged up to the wavenum- 
ber k ~ 0.3/Aei and is enhanced at higher wavenumbers (Fig. ITUb). It can reach values up to 
an order of magnitude higher than the damping rates given by the usual collisional approach 
[20I 1 close to the short wave length applicability limit of our model. In a low-/3 plasma, the 
slow MHD wave and the ion-acoustic wave can experience enhanced damping due to nonlocal 
effects (Fig. ITTb) while the fast MHD mode has the usual damping. The enhanced damping 
of the slow MHD wave appears at A; > O.l/Agj and remains up to A; ~ 3/Xei until the slow 
MHD wave becomes a second ion cyclotron wave. At the same wavenumber, k ~ 3/Aej, the 
ion-acoustic wave appears (the curve 2 in Fig. ITTj) . The damping of ion-acoustic wave can 

)e much higher than the standard collisional damping. This agrees with the result of Ref. 

3] because of the neghgible role that is played by the magnetic field in this case. 



VII. RELAXATION OF HOT SPOT TEMPERATURE IN A MAGNETIC FIELD. 

Current laser plasma interaction experiments involve random phase plate laser beams 
which control the intensity distribution in the focal plane. This makes a hot spot geometry 
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for the temperature inhomogeneity a basic element of transport theory in such plasmas. 
Typically the characteristic scales of the temperature gradient are in the regime of nonlocal 
transport. A self-generated magnetic field has also been well documented in such plasmas. 
Its existence motivates the discussion below, where we apply our transport theory to a single 
hot spot transport problem with magnetic field effects included. 

We consider the relaxation of a single hot spot in a plasma with a given magnetic field. 
The geometry of our model involves an initial temperature perturbation with cylindrical 
symmetry which evolves due to thermal transport. The magnetic field is assumed orthogonal 
to the radial temperature gradient, i.e. it can have only Bz or components. We will 
examine conditions in laser produced plasmas where magnetic field effects influance transport 
in the nonlocal regime. The solution of the initial value problem follows Ref. [22], where the 
relaxation of a hot spot was studied in an unmagnetized plasma. A temperature perturbation 
in a cylindrical hot spot is assumed to have a Gaussian profile 

dT,{r)=5Toexp(-^] (31) 



at time, t = 0. We also assume that the plasma flow has a negligible effect on the temperature 
relaxation, u = 0. 

Our model includes the energy balance equation for electrons 

d6T, . 2 



—I- 



-kx ■ q± , (32) 



dt 3ne 

where the heat flux across magnetic fleld, q_L, (fT^ has the following form ■ = 
—i\<.\K,±6Te. The value n± = x± ~ Te{a''^'^ / a ±_ is the transverse nonlocal heat conductivity. 



22| . The radial temperature 



The solution to the equation (jH^ is given by Eq. (5) in Ref. 
and electron heat flux proflles are presented in Fig. ^Jfor the magnetic fleld corresponding 
to Hall parameter fi/i/^ = 0.3. The result of the classical transport approach (dashed 
line) overestimates the electron heat flux (Fig. IT^). Because of this reason, the hot spot 



temperature obtained from local theory 16[ decreases faster (Fig. [T^ l than predictions 



from Eq. (jozj). As discussed in the previous section, a combination of the magnetic fleld 
and nonlocal effects leads to a slightly larger heat flux and smaller temperature as compared 
to the nonlocal and unmagnetized case. This is illustrated by dotted (present theory) and 
continuous lines (nonlocal and unmagnetized theory from Ref. js^) in Fig. ^| 

As expected, effects of nonlocal thermal transport become negligible with increasing 
values of the magnetic fleld. This is illustrated by Fig. ^1 where the relaxation time 
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(the time necessary for the temperature in the center of the hot spot to drop to half of 
its initial value) is shown as a function of the hot spot size, R/Xei for different magnetic 
fields. For R/Xei ^ 3 and fi/z/^ = 1 the effect of nonlocality almost vanishes. For small Hall 
parameters, ^/lyJi ^ 0.1, the influence of a magnetic field on the temperature relaxation can 
be ignored. For the intermediate magnetic field strength 0.1 < ^/i^Ji < 1 both magnetic 
field and nonlocal effects are important and their interplay defines the relaxation of the hot 
spot temperature. 



VIII. CONCLUSIONS 

We have applied a systematic closure procedure [3| to a magnetized plasma in the 
regime of moderate collisionality where fc||Aei < 1 and k± < Max{l/p,l/Xei}- We have 
derived nonlocal hydrodynamic equations which are valid for small perturbations about 
an equilibrium homogeneous state and for a plasma with high ionic charge. Our transport 
theory is fully equivalent to the linearized electron kinetic equation but possesses advantages 
of the reduced model which describes plasma in terms of just few moments of the distribution 
function. The closure procedure provides nonlocal relations between electron fluxes and 
gradients in the plasma density and electron temperature, and an electrostatic field. These 
tensor transport relations account for all orientations of the magnetic field and the plasma 
gradients: longitudinal, oblique, and transverse. 

Competition between nonlocal and magnetic field effects is poorly understood and our 
theory provides a rigorous means of examining the relative importance of both factors. We 
have found that large magnetic fields contribute to the enhancement of electron conduc- 
tivity in the weakly collisional regime. This increase can offset inhibition due to nonlocal 
effects. This is an important effect that is necessary to explain transport phenomena in laser 
produced plasmas p|. The practical importance of our theory lies in providing benchmark 
results for comparisons with the phenomenological models such as Ref. ^ or future kinetic 
simulations (cf. Refs. 13, 28^). 

We have described the application of nonlocal hydrodynamics to the MHD wave disper- 
sion and damping calculations and to hot spot temperature relaxation. There are important 
enhancements to collisional damping of MHD modes in the nonlocal transport regime. The 
presence of a strong magnetic field is the main source of thermal transport inhibition in the 
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relaxation of a hot spot. For the intermediate values, 0.1 < ri/z/^ < 1, the hot spot displays 
the simultaneous effects of electron gyro motion and nonlocal transport. 

Our results provide a next step in the ongoing efforts to develop a plasma hydrodynamics 
model which is valid for the wide range of plasma conditions and includes nonlinear re- 
sponse to large am plit ude perturbations. The application of linear theory to experimental 
data, such as 0, 0, 24 , , and comparison with kinetic simulations and its generalization 
to include local dependence on plasma parameters, as in Ref. could lead to a nonlocal 



nonlinear transport theory in magnetized plasmas. Linear nonlocal hydrodynamics models 
provide a description of the plasma thermal response to electromagnetic radiation in stimu- 



lated Brillouin scattering, filamentation and in the filament resonant instability 



23 



Current results will allow generalization of these processes to the case of magnetized plasmas. 
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APPENDIX A: NONLOCAL TRANSPORT COEFFICIENTS 

The components of the electrical conductivity tensor, ct, and thermocurrent tensor, a.^ 
are presented in the following form: 



0-11 = < - 0-± + (Xf ^ ^ , (712 = 0-° - o-^ - ' 



rr 4- 



0-21 = cr^ - (Ja + cTa't? ' ^22 = (^± - o-± + 



"11 = a" - + af'^ , ai2 = a° - - af'^ , (Al) 

"21 = "a - "a + "a'p ? a22 = a±- a]_ + a;,^— , 

g El W E2 

"31 = «|| - aJi - "ii p , «32 = -p • 
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The coefficients cr° and q;° are given by the transport theory for strong coUisions ()16p and 
the introduced new transport coefficients, such as 



a 



r)TN 



El 



J 



r)NT 



a 



E2 



erif. 



D 



a 



NT 
qaT 



Xa 



-L el^ ^ NT 

riNT r,TN 



Ha 



J^NT 
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dEiN 
^TN 


r)NT 


PfNT 
^NT 




r,E2T 
J-^NT 




r,NT 


^NT 


Jq 

Ha 


r,EiN 
J^TN 


r)NT 

^NT 


r,NT 
J^NT 



DE2N 

^TN jT 
DNT '^a 
^NT 



J 



fc2 



jjNT 
J^NT 



a 



E2 



Tek^ 



p,E2T p,E2N 
jE2 _ ^NT jN _ ^TN tT 

'^Qa riNT '^qa DNT '^qa 



(A2) 



'NT ^NT 

Expressions in ()A2|) can be calculated in terms of moments of the basic distribution functions 
satisfying Eq. (jSI). 



APPENDIX B: DISPERSION RELATION 



This appendix presents the result of calculation of the dispersion relation for MHD waves. 
We write the nonlocal hydrodynamic equations (j28j) and the wave equation (j29|) in the matrix 
form. The vanishing of the following determinant leads to the dispersion relation. 



Pa 



k^aT^ 






kl 
kl 

k±a^^Te 
-ukjc^/AiT 









ik'^a^ 











-CTa- Cri2 


ik']_a''^ 


-1 










—ik\a\ 





-1 





-k±aiTe 















i 


[^'- kW^ 









4:71 














iiu"^— k'^v\) k I VL. 



ik± 

eUe 



Ane 



k^cl 



where {a, = k '^{k'^,a\\ + k'j_a±, k?,a\\ + k1a±, k?,x\\ + kj_x±}, Cs = \/ ZTe/rrii is the 



ion-acoustic velocity, and va = B / Aumirii is the Alfven velocity. 

This dispersion relation could be rewritten in the form (jHUj) where the parameters A, B, 
C, D, and E, which are the functions of the nonlocal transport coefficients arise naturally. 
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They are 

A = ne{k\uja3 + iuj'^k'^a-) /2 , 

B = rieOj'^ik'^uja — ik\aJ) 

C = k\k'^{x(^- ~T(,a ) + iuj'^k'^{x'^ ~TeCi^)-i 

D = ujk^{xa+ - Tea+) - ik\k^{xcr3 - Tea^), 

E = —k']uj(k'^(aa/\ — ckta) + ik^cuas) , 
e 

where we introduced the following combinations of the transport coefficients: 

k^ 

a_ = (j||cr_L + aan , o"+ = cr{a± + (T22) + "p^'^A' 

kj kl 
a- = (— ay(T_L + -p-aio-y + a^fTii), 

a+ = a^((T_L + 0-22) + -p'(2a±aAC^A - "aC^), 

as = (2cr||(aAa± - (T±a/,) + 2(Tii((TAa - o-Qa) + (0-12 + ■^<^2i){(^\\a± - (^±a\\)), 

kl k'i 

ft- I / , s , s ,2 I II 2 2 ' II 



04 = -^a/s,{als + cr_) + {a± + cr22)«- + (a (Ti2(T2i + p«||'^A + ««||0"a(o'i2 + pO"2i)), 



kf. 



k^ 

The derived dispersion relation (j3Up describes small amplitude MHD modes in the k\\Xei < 1- 
range of the wave numbers. This encompasses the range of wave numbers from the adiabatic 
perturbations for small k to the isothermal perturbations for shorter wavelength. 
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Figure captions 

FIG. 1: Longitudinal electrical conductivity normalized to (Jq = 32e^ne/37rmef^ for k±/k\\ = 1000 
and n/iyj- = 0.01 (1), 0.1 (2), 0.3 (3), 1. (4), 3 (5), and 10 (6). 

FIG. 2: Longitudinal temperature conductivity normalized to xo = (200/37r)nei^Te-^ei for k±/k\\ = 
1000 and n/iyei = 0.01(1), 0.1 (2), 0.3 (3), 1. (4), 3 (5), and 10 (6). 

FIG. 3: Longitudinal electrical conductivity normalized to ctq = 32e^ne/37rmei'^ for k± = fcy and 
n/uj. = 0.01(1),0.1 (2), 0.3 (3), and 10 (4). 

FIG. 4: Longitudinal temperature conductivity normalized to xo = (200/37r)nefreAei for k± = fey 
and n/i^ei = 0.01 (1), 0.1 (2), 0.3 (3), and 10 (4). 

FIG. 5: Transversal electrical conductivity normalized to ctq = 32e^rae/37rmet'2 for k± = k\\ and 
O/j/J = 0.01(1),0.1 (2), 0.3 (3), and 1 (4). 
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FIG. 6: Transversal temperature conductivity normalized to xo = (200/37r)neUTeAej for k± = fey 
and n/i^ei = 0.01 (1), 0.1 (2), 0.3 (3), 1 (4) and 3 (5). 

FIG. 7: Oblique conductivity normalized to e^ne/niei'^- for kj_ = k\\ and ^/v'^i = 0.01(1), 0.1 (2), 
0.3 (3), 1. (4), 3 (5), and 10 (6). 

FIG. 8: Oblique temperature conductivity normalized to UeVTe^ei for k± = ^|| and ^/i^'^i = 0.01(1), 
0.1 (2), 0.3 (3), 1. (4), 3 (5), and 10 (6). 

FIG. 9: The ratio of the transversal temperature conductivity to the longitudinal temperature 
conductivity for kXd = 0.1 (large dots) and 0.3 (small dots). Solid line corresponds to the local 
Braginskii's theory. 

FIG. 10: Solution (dots) to the dispersion equation for the MHD waves for real (w) and imaginary 
(7) parts of the frequency (panels a and b, correspondingly) in comparison with the result of 
conventional theory (solid lines) for a plasma with p = 50 and h = 10. 

FIG. 11: Solution (dots) to the dispersion equation for the MHD waves for real (w) and imaginary 
(7) parts of the frequency (panels o and b, correspondingly) in comparison with the result of 
conventional theory (solid lines) for a plasma with p = A and /i = 50 

FIG. 12: The temperature (a) and heat flux (6) profiles (dots) at the time t = 2/i^^ in comparison 
with the results of nonlocal theory without magnetic field (solid lines) and classical local 
approach Q| (dashed lines). The plasma parameters are: ^/i^'^i = 0.3, R/\ei = 3, and Z = 5. 

FIG. 13: The temperature relaxation time as a function of the hot spot radius (dots) in comparison 



with the nonlocal theory without magnetic field 
[1^ (dashed lines) for different Hall parameters. 



221 (solid lines) and the classical local approach 
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